Real-time determination of free energy and losses in optical absorbing media 
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We introduce notions of free energy and loss in linear, absorbing dielectric media which are relevant 
to the regime in which the macroscopic Maxwell equations are themselves relevant. As such we solve 
a problem eluded to by Landau and Lifshitz [|l|L in 1958, and later considered explicitly by Barash 
and Ginzburg Q, and Oughtsun and ShermaiiQ. As such we provide physically-relevant real-time 
notions of "energy" and "loss" in all analogous linear dissipative systems. 
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In a previous publication we showed that in many in- 
stances fast and slow light ||, |6| is a manifestation of 
a dielectric interacting differently with the early parts of 
an electromagnetic pulse than with its later parts . For 
example, slow light in a passive linear dielectric typically 
corresponds to energy in the leading part of the pulse be- 
ing preferentially stored by the medium and then being 
largely returned to the pulse's (backward) tail. Quanti- 
fying the extent to which this process is possible was a 
primary motivation for this work. 

Application of the principles presented in this letter 
will in no-wise supersede the group velocity description 
of such phenomena. The present development compli- 
ments other such analysis by attaching precise notions of 
reversibility and irreversibility to the medium's storage of 
the pulse energy. As such it also solves a long outstand- 
ing problem regarding the real-time meaning of " energy" 
and "loss" in linear dissipative systems. 

From Foynting's energy conservation theorem (Eq. |^) 
the total energy density in a dielectric at time t, u{t), is 
the sum (Eq. ||) of the field energy Ufieid{t) (Eq. 0) and 
the medium- field interaction energy Uint.^t) (Eq- m- 



u{t) = Uficld{t) + Uint.{t) 



Ufieldit) 



Utnt.it) 



E{t)P{t) dr. 
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Here we restrict to isotropic, temporally dispersive me- 
dia in which a scalar analysis suffices, and in which one 
may safely suppress reference to the spatial coordinate 
X. In addition we restrict to linear media. Consequently 
we restrict to the case in which the medium's response to 
an applied field E is completely determined by a scalar. 



point-wise defined susceptibility x — x(ti',x). For a ba- 
sic discussion of properties of the medium-field interac- 
tion energy Uint.it) anisotropic media, see 1^. Since 
Uint.it) quantifies the net work the field has done against 
the polarized medium in the course of creating the cur- 
rent system state, in the following we will also call Umf . (t) 
the (medium) internal energy. In the sequel, notions of 
work done, and the ability to do work in the future, will 
become pivotal in establishing an unambiguous, dynam- 
ically relevant notion of energy allocation in dielectric 
media. 

Landau and Lifshitz interpreted Mi„(.(-|-oo) as the en- 
ergy that is eventually lost to and dissipated by the 
medium over the course of the medium-field interac- 
tion [Q. This asymptotic quantity only depends on the 
medium susceptibility x ^-nd the evolution of the elec- 
tric field E, not on the particulars of any microscopic 
model giving rise to x- This suggests the possibility 
of establishing a real-time, model independent notion of 
loss. Barash and Ginsburg considered this question 
They concluded, however, that a real-time determina- 
tion of losses is impossible without a microscopic model, 
i.e. that a knowledge of the macroscopic susceptibility 
X is insufficient. Having concluded that a model inde- 
pendent notion of loss is meaningless, they calculated (a 
certain notion of) loss for specific models. In particular 
they generalized the work of Loudon |^ concerning sin- 
gle Lorentz oscillator media by making a straightforward 
extension of his notion of losses to multiple Lorentz oscil- 
lators. Loudon's and, so, Barash and Ginsburg's notion 
of losses amounts to identifying those terms in (a certain 
expression of) the internal energy which explicitly de- 
pend on phenomenological damping parameters. More 
convincingly, their determination of the "energy" in the 
medium amounts to summing the kinetic and potential 
energies of the individual oscillators. We call this notion 
of energy the polarization energy, Upoiarizationit) ■ 

We initially followed a certain extension of this rea- 
soning. We hoped to establish a unique map from an 
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FIG. 2: The instantaneous losses of a two oscillator (long 
dash), and a coupled oscillator (short dash) model of a double 
Lorentz oscillator susceptibility responding to an instigating 
delta function E-field at time t = 0. The interaction energy 
'u.int\E\{t') (solid), and the irreversible energy Uirrev.{t) (dot- 
dash) are also shown. 



FIG. 1: Macroscopic equivalent spring-mass systems for a 
double Lorentz oscillator and their LRC circuit analogs, (a) 
The two oscillator representation, (b) The LRC circuit analog 
of (a), (c) The coupled oscillator representation, (d) The 
circuit analog of (c). 



arbitrary susceptibility xi'-^) to an oscillator representa- 
tion, which in turn would establish a unique (polariza- 
tion) energy, as well as losses via Loudon's identification 
procedure. Unfortunately, we found that generically a 
susceptibility can be mapped to many different micro- 
scopic models. This property is perfectly analogous to 
the one by which it is possible for distinct LRC circuits 
to yield identical impedances. However, we then con- 
jectured that all such equivalent representations might 
yield identical values for the polarization energy. Unsur- 
prisingly, one finds that this is not the case dynamically. 
(Although asymptotically all such representations must 
obviously agree: Upoiarizationi+oo) = 0, and the losses 
are given, then, by Ui„t{+oo), which, as mentioned, is 
completely determined by %.) 

Figures ^ and ^ demonstrate this ambiguity. Given 
a double Lorentz oscillator susceptibility x{uj) ~ 
Sn=i 2^p„/(^n ~ ^Tn"^ ~ ^'^) examine two different 
microscopic models for the same susceptibility. The first 
model Xa(w), is given by the explicit structure of x{^) 
as just written. It corresponds to two different masses of 
equal charge on two different springs, each spring hav- 
ing different restoring and damping parameters. The 
total polarization is given by the net displacement of 
both charged masses from their equilibrium positions: 
P = Xi^+ X2„ = XaE. See FIG. |l|(a). 

We use the tangency algorithm jl^, a transfer func- 
tion preserving algorithm intended for model reduction in 
LRC circuits, to construct the second microscopic model 



Xb{^)- Mechanically this model corresponds to coupled 
oscillators. See FIG. |](c). The susceptibility for this 
microscopic model is given in Eq. (^ (where we have 
mapped to to is for convenience). In this model the total 
polarization is given by the displacement of the charged 
mass from its equilibrium position P = Xi^ = XbE. It 
can be shown that the coefficients fc„^ , 7„^ , can be 
found such that Xai'^) = Xb{^) for all frequencies to. 



By inserting the two representations into the internal 
energy Uint{t) = E{t)P{t) dr one can find the ki- 
netic and potential energies associated with each rep- 
resentation. This polarization energy, the energy reck- 
oned instantaneously via the motion and position of the 
masses, will be shown to be representation dependent. 
To demonstrate this, we calculate the "microscopic" en- 
ergetics associated with the response of these macroscop- 
ically equivalent systems to an impulse. That is we excite 
the two media, with susceptibility models Xa and Xb, via 
a delta function E-field at time t — 0, and then plot the 
losses for each representation. We assign values to the 
parameters of the two oscillator, Xa representation (the 
plasma frequencies, resonant frequencies, and damping 
coefficients) , and then find parameters for the coupled os- 
cillators', Xb representation such that the representations 
are macroscopically equivalent. Figure ^ shows the evolu- 
tion of the losses for the different representations, these 
losses determined in the sense of Loudon, Barash and 
Ginzburg. The corresponding differences between the in- 
ternal energy imparted to the media by the delta func- 
tion E-field (as shown by the piecewise constant curve 
in figure ||) and the plotted losses gives the polarization 
energies for the two representations. 

The qualitative features of Figure ^ can be understood 
intuitively: the delta excitation of the medium instan- 
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taneously creates the polarization energies, which, then, 
decrease as the systems dissipate these energies. How- 
ever, as clearly seen in the figures, the polarization en- 
ergies for each microscopic representation differ signifi- 
cantly. Thus the polarization energy does in fact depend 
on the microscopic model ascribed to the susceptibility 
x(tj). Of course, since Barash and Ginzburg specified no 
underlying macroscopic physical principle in their deter- 
mination of "energy", this result is just a consequence 
of their (lack of a) definition of energy. Consequently, 
we will argue that a model-dependent (e.g. polarization) 
energy is irrelevant, both from the point of view of the 
relevant physical principles that should be required of a 
viable macroscopic description, as well as from that of the 
practical considerations of the energy storage and return 
process mentioned at the beginning of this article. 

Before establishing relevant macroscopic principles of 
energy allocation, we highlight the connection between 
mechanical and electrical oscillators and the ambiguities 
associated with representation in the context of the lat- 
ter. The perfect analogy between electrical LRC circuits 
and linear passive dielectric media (where polarization 
P, susceptibility x and E-field E , translate to charge Q, 
derivative of admittance dtA, and electromotive force £, 
respectively), allows one to reinterpret the polarization 
energy evolutions implied in Fig. |2| as the evolutions of 
the energies contained in the dispersive elements (induc- 
tors and capacitors) of admittance-equivalent circuits. In 
the case of electrical circuits, the plotted losses corre- 
spond to the losses in the dissipative elements of the two 
circuits, i.e. the losses through the resistors. This time 
we interpret the discrepancies in these curves as indicat- 
ing that different LRC circuits with the same admittance 
can, temporarily, lose energy to their resistors at different 
rates. Figures ||(b) and (d) give LRC circuits correspond- 
ing to the mechanical models shown in Figures |l|(a) and 
(c), respectively. 

From the discussion above, one concludes that the po- 
larization energy and the associated losses depend intrin- 
sically upon the specific microscopic model giving rise to 
the susceptibility x- However, these allocations cannot 
be relevant macroscopically: for example, within the phe- 
nomenological framework in which x is introduced, the 
spatial and temporal evolution of the various fields de- 
pend only upon this particular piece of information, not 
upon its various representations. In this circumstance 
in which the system dynamics is completely determined 
by some piece of information (e.g. x)j to say that some 
other piece of information is important in order to es- 
tablish some particular notion of energy allocation is to 
admit that that notion of energy allocation is irrelevant 
to those all-ready specified dynamics. We establish a no- 
tion of energy allocation that is determined uniquely by 
X and, so, is relevant to system dynamics. In particular 
it provides a precise notion of the maximum possibility 
for the field to recover energy from the medium in, for 



example, the borrow-return process mentioned at the be- 
ginning of this article. 

In this article, we introduce the irreversible energy den- 
sity Uirrev.{t)- At any given time it is defined to be the 
infemum of all possible [LL] asymptotic losses Mi„t.(-|-oo), 
this extrema being realized over all possible future evo- 
lutions of the electric field i?, holding its past evolution 
fixed. Temporarily using notation emphasizing that the 
internal energy at time t depends not only on this time, 
but also on the history of the instigating electric field E 
up until that time, one writes 

Uirrev.[E]{t) = Uirrev .[E^^]{t) := iuf + E^]{+Oo). 

(6) 

Here E^ denotes the electric field time series E{t) with 
its i-future (t > t) eliminated. Similarly, E^ denotes an 
appended electric field time series E{t) with its i-past 
(r < t) eliminated. For a passive dielectric, this infemum 
exists and, so, is unique [0. In particular, and is shown 
later in Eq. (^, it does not depend upon an explicit 
representation for x- 

Almost tautologically, from definition (||), it follows 
that Uirrev. (t) Can never decrease as time increases. Thus, 
at any given time t, it quantifies a component of the 
medium internal energy Uint.(t) that will, under all cir- 
cumstances, remain in and be dissipated by the medium. 
That is it specifies a component of the medium inter- 
nal energy that cannot under any circumstances be re- 
turned to the field. Moreover, since it is defined in terms 
of an infemum, i.e. a greatest lower bound, all notions 
of loss greater than this value are too pessimistic: at 
any given time t, there always exists a future medium- 
field interaction creating less eventual energy loss to the 
medium than any value greater than that specified by 
Uirrev. {t) ■ This is truc regardless of how small this overes- 
timation is. Consequently, within the phenomenological, 
macroscopic framework in which x dictates the system 
dynamics, this quantity uniquely records the energetic 
irreversibility generated within this dissipative system. 

In Fig. ||the dot-dashed curve specifies the irreversible 
energy for the case considered, obviously valid for ei- 
ther of the two x-equivalent physical systems. Note 
that in the figure it is never exceeded by the losses in 
the sense of Loudon, Barash and Ginzburg. Indeed one 
can prove that this relationship must hold between the 
macroscopically relevant irreversible energy and any no- 
tion of loss specified microscopically. Equivalently one 
determines that energy (the ability to do work) speci- 
fied macroscopically cannot exceed any such microscopic 
notion. Further, one concludes that the former is almost 
always strictly less than the latter because of incoherence 
among the system's microscopic, energy containing ele- 
ments. The former statement (the one regarding losses) 
can be obtained by repeated application of the following 
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theorem: (As in the monotonicity of Uirrev.{t)^ theo- 
rem follows almost tautologically from the definition.) 

[E;xi+X2m > [E;xim + [E;x2]{th 

(7) 

with strict inequality holding almost always for non- 
trivial xi Etnd X2- Eq. (^) demonstrates that in "ad- 
ditive" processes (i.e. generating media mixtures), ir- 
reversibility is generated. By repeated (additive) sub- 
divisions of a macroscopically relevant susceptibility x 
into pieces Xij one may obtain microscopic representa- 
tions of the medium response. If the elements are "sim- 
ple" enough (to be quantified later) Uirrev.[E;xi] will be 
equivalent to the parameter-dependent notion of loss sug- 
gested by Loudon, and Ginzburg and Barash. Indeed the 
notion of loss and energy specified by Loudon happen to 
agree with the macroscopically/phenomenologically rele- 
vant notion herein introduced in the cases he considered- 
single Lorentz oscillator media. [This is not the case 
for (competing) multiple oscillator media considered by 
Barash and Ginzburg, as demonstrated by figure |^.] 

The difference between the current internal energy and 
the current irreversible energy gives the reversible energy 

^rev. (^) ■ 



.it) 



it) 



.it)- 



(8) 



The reversible energy gives the least upper bound on 
the amount of energy that the medium can relinquish af- 
ter time t: any amount greater than this value, no mat- 
ter how small the discrepancy, overestimates the ability 
of the medium's microscopic, energy-containing compo- 
nents, say, to organize themselves and do useful, macro- 
scopic work. Obviously the difference between the (piece- 
wise constant) internal energy plotted in Fig. ^, and the 
dot-dashed curve in that figure, gives this dynamical no- 
tion of the possibility for the medium to do work (against 
the field) for the example considered. 

Using Eq. (||), one immediately shows that Uint.it) 
is constant after the electric field ceases. From defini- 
tion (^), and the theorem regarding the monotonicity 
of Uirrev.it), it follows that Urev.it) cau ucver increase 
after such time, i.e. after the electric field quits subsi- 
dizing its existence by doing work against the polarized 
medium. We will show that Urev.it) is never negative. 
Consequently, one sees that when the system becomes 
energetically closed, the reversible energy behaves like a 
dynamical system free energy (density) , equivalently like 
a system Lyapunov function (density). For this reason, 
and for the microscopic consideration mentioned, in par- 
ticular the entropy-generation-like statement embodied 
in Eq. (^) , we will also designate the reversible energy as 
the (medium-field) free energy. 

We finish with a formula demonstrating how the 
macroscopic loss (and so the free energy (via 0)), can 



be calculated. (In particular how the irreversible energy 
plotted in figure (||) can be generated.) This formula is 
obtained by applying a variational principle to the defini- 
tion (^ , and solving the resulting Riemann-Hilbert prob- 
lem. One finds that, for passive, causal dielectrics. 
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where 



A — lim 
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(10) 



Here we introduce the instantaneous spectrum at time 
t, Etiuj). It is the Fourier transform of (see the lines 
following definition (^). We also introduce the medium 
complexity factor (/)-|_(w). Its deviation from unity gives 
a measure of the effective macroscopic/phenomenological 
incoherence of possible microscopic, energy containing el- 
ements. Media for which (f>+iLij) is identically unity we 
call simple media, the rest we call complex. In the case 
that the susceptibility x corresponds to a single Lorentz 
oscillator medium, as in the case considered by Loudon, 
4>+iuj) is identically one, and Eq. (I) reduces to 



u.[E]{t)^\ / P^t) dr 



(11) 



where, then, A is determined in terms of the phenomeno- 
logical damping parameter 7 and the effective plasma 
frequency LOp-. 



A = 



(12) 



In such case, then, and as claimed by Loudon, the rel- 
evant losses correspond to the frictional losses generated 
by the single Lorentz oscillator. 



SUMMARY 

We have introduced notions of free energy and losses 
relevant to the macroscopic behavior of passive, linear 
dielectric media. These notions are relevant in the same 
regime in which the macroscopic Maxwell equations are 
themselves relevant, i.e. in the regime in which that the- 
ory specifies all measurable dynamics. In particular, the 
macroscopic theory introduced is relevant to the regimes 
in which the macroscopic, energy barrow-return process 
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describes the production of slow and fast light in pas- 
sive, linear, temporally dispersive media. Further com- 
munications will describe the precise evolutions of the 
ideal medium-field interactions giving rise to maximum 
energy recovery, i.e. those "recovery" fields suggested by 
the variational definition . The nature of the analogy 
of such recovery fields in dissipative media with reversible 
processes in the thermodynamic setting will be analyzed. 
Finally, the nonintuitive evolution of the medium- field in- 
teraction on such recovery fields in complex media will 
be exposed, e.g. the failure of monotonicity in the evo- 
lution of the internal energy to its minimum value, even 
on ideal recovery fields, and the identification of such as 
a measure of the level of macroscopic disorder created by 
the existence of many, competing degrees of freedom. 

Conversations with Kurt Oughstun, Joseph Eberly, Pe- 
ter Milonni, Aephraim Steinberg, Raymond Chiao, and 
Michael Fleischhauer are gratefully acknowledged. This 
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Foundation under Grant No. PHY99-07949. 



pp. 274-332. 

[2] Yu S. Barash and V. L. Ginzburg, Usp. Fiz. Nauk 118, 
523-530 (1976). [English trans.: Sov. Phys.-Usp. 19,263- 
270 (1976)]. 

[3] K. E. Oughstun and G. C. Sherman, Electromag- 
netic Pulse Propagation in Causal Dielectrics, (Springer- 
Verlag, Berlin, 1994), pp. 14-31. 

[4] C. G. B. Garrett and D. E. McCumber, Phys. Rev. A 1, 
305-313 (1970). 

[5] S. Chu and S. Wong, Phys. Rev. Lett. 48, 738-741 (1982). 

[6] L. P. Hau, S. E. Maris, Z. Button, and C. H. Behroozi, 
Nature 397, 594-598 (1999). 

[7] J. Peatross, M. Ware, and S. Glasgow, J. Opt. Soc. Am. 
A 18, 1719-1725 (2001). 

[8] S. Glasgow, M. Ware, and J. Peatross, Phys. Rev. E 64, 
046610 (2001). 

[9] R. Loudon, J. Phys. A. 3, 233-245 (1970). 
[10] S. Glasgow, M. Meilstrup, M. Clayton, G. Hovhannisyan, 

and C. Broadbent, (unpublished). 
[11] C. Broadbent, G. Hovhannisyan, M. Clayton, J. 
Peatross, and S. Glasgow, AMEREM 2002 Conference 
Proceedings, (accepted for publication). 



[1] L. D. Landau and E. M. Lifshitz, Electrodynamics of 
Continuous Media, 2nd ed. (Pergamon, New York, 1984), 



